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We investigate the effect that the temperature dependence of the crystal structure of a two 
dimensional organic charge-transfer salt has on the low-energy Hamiltonian representation of the 
electronic structure. For that, we determine the crystal structure of k-(BEDT-TTF)2Cu2(CN)3 
for a series of temperatures between T = 5 K and 300 K by single crystal X-ray diffraction and 
analyze the evolution of the electronic structure with temperature by using density functional theory 
and tight binding methods. We find a considerable temperature dependence of the corresponding 
triangular lattice Hubbard Hamiltonian parameters. We conclude that even in the absence of change 
of symmetry, the temperature dependence of quantities like frustration and interaction strength can 
be significant and should be taken into account. 

PACS numbers: 71.15.Mb, 75.10.Jm, 61.05.cp, 61.66.Hq 



I. INTRODUCTION 

The two-dimensional organic charge-transfer salts 
based on bis(ethylenedithio)tetrathiafulvalene (BEDT- 
TTF or even shorter ET) molecules in a K-type lat- 
tice arrangement have been intensively studied over the 
past thirty years due to their complex interplay between 
electron correlation and the effects of low dimensional- 
ity and spin frustratiouii In particular, the discovery 
of spin-hquid behavior in K-(BEDT-TTF)2Cu2(CN)3l 
has fascinated experimentalists and theorists alike. Is- 
sues of current interest concern the nature of the low- 
temperature spin-liquid realized in this materials."'^, and 
the various anomalies observed upon approaching the 
spin-liquid state from high temperatures. These anoma- 
lies include drastic changes in the -"^H-NMR relaxation 
rate around 200-150 K— , the thermopower at 150 K^, 
relaxor-type ferroelectricity around 60 and a mysteri- 
ous phase-transition anomaly at 6 . The latter feature, 
which manifests itself in anomalies in thermodynamio^i^ 
and transport^ quantities, is accompanied by pronounced 
lattice effects.^ Various scenarios have been suggested for 
the 6 K anomaly including a crossover from a thermally to 
a quantum disordered stated, an instability of the quan- 
tum spin-liquid^' or a distinct type of charge order- 
ingjil Theoretically, the spin-liquid properties have been 
investigated on the basis of the anisotropic triangular- 
lattice Hubbard HamiltonianJ^ii^r— The parameters t, 
t' and U of this Hamiltonian have been determined with 
semi-empirical^ as well as first principles methods^^'^^ 
based on the experimental structure at room tempera- 
ture. Missing in such investigations is, however, the con- 
sideration of a possible temperature dependence of the 
model parameters. 

In this work, we find by a combination of single crys- 
tal X-ray diffraction at various temperatures and density 



functional theory calculations that even in the absence 
of structural phase transitions, the temperature depen- 
dence of the structural parameters is significant enough 
to influence the electronic behavior and the determina- 
tion of the Hamiltonian model parameters in k-(BEDT- 
TTF)2Cu2(CN)3. We suggest that this has subtle effects 
on the degree of frustration and interaction strength. 



In 1991, the crystal structure of k-(BEDT- 
TTF)2Cu2(CN)3 was first reported^ and the unit 
cell parameters confirmed on twinned crystals at 
room temperature. The crystal structure was rede- 
termined in 1993, also on a twinned crystal at room 
temperature;^ The crystal structure of the very similar 
k'-(BEDT-TTF)2Cu2(CN)3 structure, which is reported 
to be an ambient pressure superconductor, was described 
in 1992<22iS^ The crystal structure of the k phase was 
subsequently redetermined in 1997,^^ and the unit cell 
redetermined in 2001, both at room temperaturei^ 
Already in 1991, the unit cell parameters were reported 
at 300 and 30 K, with a(^30K)/a(300K) = 0.9964, 
b{30K)/b(300K) = 0.9932, and C(^30 k)/c(^300 k) = 0.9900, 
and /3(30if)//3(300K) = 1.0146-^21 Our compara- 
ble contraction values between 300 and 20 K are 
a{20K)/ 0,(300 K) = 0.9988, &(20if)/^(300K) = 0.9958, and 
C{20K)/c(3ooK) = 0.9905, and 13 (20 k)/I3 {300 k) = 1.0148 
and agree well with the previous results. In addition, 
we find evidence for an ordering of the ethylene groups 
in a staggered conformation between 200-150 K. Despite 
the current interest in this material as a spin-liquid 
candidate, no low-temperature structural determinations 
have yet been reported. Herein, we present a detailed 
characterization of the crystal structure as function 
of temperature determined on a single crystal with 
agreement factors better than any previously reported. 
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II. SINGLE CRYSTAL X-RAY DIFFRACTION 

A black, plate-like crystal of k-(BEDT- 
TTF)2Cu2(CN)3 with dimension 0.4 x 0.4 x 0.4 mm^ 
was placed onto the tip of a glass fiber and mounted 
on a Bruker APEX II 3-circle diffractometer equipped 
with an APEX II detector. Temperature control in the 
100-300 K region was provided by an Oxford Cryostream 
700 Plus Cooler while below 100 K it was provided 
by a Cryocool-LHE cryogenic system (Cryo Industries 
of America). The sample temperature below 100 K 
was confirmed by installing a Cernox thermometer 
(Lakeshore) in the immediate vicinity of the crystal and 
stabilizing the temperature immediately prior to data 
collection. In order to reduce adverse thermal effects the 
thermometer wires (twisted manganin wires, 0.5 mm 
diameter, supplied by Lakeshore) were anchored on the 
surfaces of the cryostat exposed to *He stream. The data 
were collected using MoKa radiation (A = 0.71073 A) 
with a detector distance of 50 mm. Unit cell parameters 
were determined upon cooling in 10 K increments 
between 100 and 290 K, with a frame exposure time 
of 10 seconds. Full data sets for structural analysis 
were collected at temperatures of 5, 20, 100, 150, 200, 
250 and 300 K. The uncertainties in the temperature 
determination are typically ±0.2 K down to 100 K, 
while for the 20 K and 5 K data points an error bar 
of ±1 K has to be accepted. These data collections 
nominally covered over a hemisphere of reciprocal space 
by a combination of three sets of exposures. Data to a 
resolution of 0.68 A were considered in the reduction. 
The raw intensity data were corrected for absorption 
(SADABS'^^). The structure was solved and refined 
using SHELXTLi^ A direct-method solution was calcu- 
lated, which provided most of atomic positions from the 
E-map. Full-matrix least squares / difference Fourier 
cycles were performed, which located the remaining 
atoms. All non-hydrogen atoms were refined with 
anisotropic displacement parameters. The hydrogen 
atoms were placed in ideal positions and refined as riding 
atoms with relative isotropic displacement parameters. 

As a representative example, in Fig. [T] we show the 
structure of k-(BEDT-TTF)2Cu2(CN)3 at a tempera- 
ture of r = 5 K. In the be plane, the BEDT-TTF 
molecules form their typical K-type arrangement. One 
of the cyanide (CN~) groups resides on an inversion cen- 
ter, thus requiring it to be disordered with a 50% carbon 
and 50% nitrogen distribution on these two atomic posi- 
tions^*. Analysis of the low-temperature structural data 
(see Tab.|T]) indicates that overall, a, b and c-axes decrease 
with temperature while the /3 angle increases monotoni- 
cally upon cooling. We find that the orthogonal projec- 
tion of the a-axis, a± — asin/3, has the greatest relative 
contraction with temperature. One ethylene group is dis- 
ordered at room temperature with a staggered conforma- 
tion 77% of the time. As the temperature is lowered, the 
ethylene group remains partially disordered down to a 
temperature of 200 K and is fully ordered in a staggered 



conformation at 150 K. The anion layer becomes slightly 
more buckled at low temperature: at room temperature, 
the Cul atom lies 0.050 A out of the C11-C12-N12 plane, 
increasing to 0.072 A at low temperature. For the pur- 
pose of performing density functional theory calculations, 
the original symmetry of P2i/c has been lowered to P2i 
by choosing one of two possible orientations of the CN~ 
group in the inversion center. 




FIG. 1: Structure of k-(BEDT-TTF)2Cu2(CN)3 at T = 5 K. 
Note that the disordered CN~ group in the inversion center 
is only shown in one conformation. 



III. TEMPERATURE-DEPENDENT 
STRUCTURAL PARAMETERS 



In Fig. m we show the evolution of the lattice pa- 
rameters with temperature. Over the large investigated 
temperature range from T = 300 K down to T = 5 K, 
the volume is monotonously decreasing with temperature 
(see Fig. [5] (a)). The monoclinic angle /3 between a and 
c lattice vectors (see Fig. [2] (b)) first rapidly increases 
upon cooling down to a temperature of T = 200 K, be- 
low which it increases more gradually. In Fig. ^ (c), the 
lattice parameters are displayed as symbols, normalized 
by their values at T = 20 K (see Tab.|T]). We also include 
the relative b and c lattice constants along the two prin- 
cipal axes as obtained by thermal expansion measure- 
ments in Refs. and [sj- The out-of-plane expansivity 
data shown there were taken in a direction perpendic- 
ular to the be plane. They are show together with the 
corresponding quantity from the X-ray diffraction mea- 
surement, a± = asin/3. 



3 





300 K 


250 K 


200 K 


150 K 


100 K 


20 K 


5 K 


a (A) 16.0919(3) 16.0848(3) 16.0781(3) 16.0703(3) 16.0746(6) 


16.072(4) 


16.062(3) 


h{k) 


8.5722(2) 


8.5749(1) 


8.5737(1) 


8.5664(2) 


8.5593(3) 


8.536(2) 


8.544(2) 


c (A) 13.3889(2) 13.3373(2) 13.2964(2) 13.2698(3) 13.2678(5) 


13.262(3) 


13.271(2) 


(°) 


113.406(1) 113.853(1) 114.273(1) 


114.609(1) 


14.852(1) 


115.088(3) 115.093(2) 


V (K-") 


1694.93(6) 1682.43(4) 1670.86(4) 1660.72(6) 1656.51(1) 






n { cf / fTin 

P (,6/'-"^ ) 


1.909 


1.923 


1.937 


1.949 


1.954 


1.964 


1.962 


fj, (mm^^) 


2.266 


2.283 


2.299 


2.313 


2.319 


2.331 


2.329 


GoF 


1.057 


1.053 


1.091 


1.157 


1.293 


1.261 


1.109 


R 


0.0311 


0.0269 


0.0238 


0.0225 


0.0220 


0.0182 


0.0204 


Rm 


0.0838 


0.0735 


0.0654 


0.0627 


0.0629 


0.0503 


0.0515 


staggered (%) 


77 


86 


93 


100 


100 


100 


100 


dCu-NNC (A) 


0.050 


0.056 


0.064 


0.069 


0.073 


0.072 


0.072 


ET tilt angle i? 


66.55 


66.25 


66.00 


65.79 


65.64 


65.53 


65.52 


C^intradimcr 


3.558 


3.538 


3.518 


3.500 


3.488 


3.470 


3.473 



TABLE I: Crystal data and structure refinement of «:-(BEDT-TTF)2Cu2(CN)3. Formula = C23H16CU2N3S16; formula weight 
Mw ~ 974.43, monoclinic, wavelength A = 0.71073 A, effective number of electrons in the crystal unit cell contributing to 
F(OOO) — 978, space group P2i/c, Z — 2. Residual factor for the reflections Ri — \\Fo\ — |-Fc||/5I1 \Po\; weighted residual 
factors WR2 = \^uj{F'^ — F'^Y / Yw{F'^Y'\''^^ \ [I > 2a{I)\\ least-squares goodness-of-fit parameter GoF = [Yl'^iFo ~ 

F'^Y /{Nd — jVp)] . Staggered (%) is the percentage of ethylene groups that are in the staggered conformation. dcu-NNC (A) 
is the distance that the Cu(I) ion is out of the plane defined by the coordinated N, N and C atoms. The ET tilt angle i? is 
measured against the be plane, dintradimcr is the orthogonal distance between BEDT-TTF dimers. 



IV. TEMPERATURE-DEPENDENT 
ELECTRONIC STRUCTURE 

In the following, we determine the electronic properties 
for the resolved crystal structures at different tempera- 
tures by employing the all-electron full-potential local or- 
bitals (FPLO)— basis. We perform all calculations on a 
(6x6x6) k mesh with a generalized gradient approxima- 
tion functionali^ In Fig. [3] we present the electronic band 
structures for the various crystal structures. In the cal- 
culation, we used the staggered (majority) conformation 
of the BEDT-TTF molecules at aU temperatures. The 
changes as function of temperature for the two bands at 
the Fermi level, corresponding to the antibonding com- 
binations of the BEDT-TTF highest occupied molecular 
orbital (HOMO) levels, are relatively small. On the other 
hand, the occupied bands down to —0.7 eV below the 
Fermi level show a significant dependence on tempera- 
ture. These bands derive from the bonding combination 
of BEDT-TTF HOMO levels and from the [Cu2(CN)3]" 
anion layer. Overall, these bands show a bandwidth that 
decreases with increasing temperature. This can be ex- 
plained by the volume increase as function of tempera- 
ture (see Fig. [2] (a)). 

Further analysis of the electronic structure requires 
a reliable identification of the bands deriving from the 
BEDT-TTF molecules. For that purpose, band weights 
have been calculated for all structures and added up for 
all atoms corresponding to the BEDT-TTF cation lay- 
ers and to the [Cu2(CN)3] anion layers, respectively. 
In Fig. m blue circles and orange triangles stand for a 



predominance of BEDT-TTF weight and [Cu2(CN)3] 
weight, respectively. This identification allows us to fit 
the BEDT-TTF derived bands to a tight binding Hamil- 
tonian 

i/TB = E^-j(4s<x+H.C.), (1) 

where {cia) create (annihilate) an electron with spin 
a at site i; the sites correspond to the positions of the 
BEDT-TTF molecules, shown as balls in Fig. El While 
a good overall fit of the band structure can be achieved 
by including six neighbor BEDT-TTF molecule distances 
up to d = 9.4 A, a very good fit describing also the 
small dispersion along the B — T direction, as shown by 
lines in Fig. |31 requires fourteen neighbor distances up to 
d = 14.5 A. 

We now proceed to analyze the temperature depen- 
dence of the tight binding parameters corresponding to 
the network of BEDT-TTF dimers that arc highlighted 
in Fig.O These dimers form a triangular lattice, with the 
hopping parameter t' connecting dimers to chains along 
the h direction and the hopping parameter t forming the 
2D connections in c direction. We are interested in esti- 
mating the parameters of the Hubbard Hamiltonian for 
the anisotropic triangular lattice 
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FIG. 2: Structural parameters of «:-(BEDT-TTF)2Cu2(CN)3 
between T = 5 K and T = 300 K. (a) and (b) show the 
volume and monoclinic angle, respectively. In (c) relative 
lattice parameters are given with the T — 20 K structure 
as reference. a± = asin/3. Symbols refer to the new data 
from X-ray diffraction (this work) while lines are measured 
thermal expansion data from Ref. 9. The experimental error 
bar is comparable to the size of the symbols. 



where < ij > and [ij] indicate summations over nearest 
and next-nearest neighbors, respectively, t and t' can be 
obtained from the molecular overlap integrals t2 to by 
considering the geometrical formulas 



(3) 



For the definition of the overlap integrals t„ — ti-j, see 
Ref. [l^. Note that inclusion of longer-range hopping into 
these formulas (for example including ^5 into t') has no 
influence on the results reported in the following. 




FIG. 3: Band structures of k-(BEDT-TTF)2Cu2(CN)3 from 
T = 5 K to T = 300 K. 
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FIG. 4: Band structure of «:-(BEDT-TTF)2Cu2(CN)3 at 
T = 20 K. Blue circles (orange triangles) indicate bands with 
a majority of BEDT-TTF and [Cu2(CN)3]~ character, re- 
spectively. The TB fit is shown with lines. 



V. DISCUSSION 

Fig. [B] summarizes our findings from the tight bind- 
ing analysis. The nearest-neighbor hopping parameters 
t, forming a square lattice, increase upon cooling down 
to a temperature of T 200 K, then decrease again 
(see Fig. H] (a)). The frustrating hopping parameters t' 
show the opposite behavior as a function of tempera- 
ture, decreasing upon cooling down to a temperature of 
T = 150 K, then increasing. Interestingly, these two ef- 
fects enhance each other when we consider their ratio t' /t 
as shown in Fig. |6] (b)). t' /t which quantifies the degree 
of frustration in the system decreases from t' /t = 0.82 at 
T = 300 K to t'/t = 0.80 at T = 150 K, then increases to 
a maximum value of t'/t = 0.86 at T = 5 K. A rough esti- 
mate for the Coulomb interaction strength U can also be 
extracted from the dimer approximation U ~ 2ti where 
ti is the BEDT-TTF intradimer hopping integral (see 
Fig. ini (c)). We find that the measure of the interaction 
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FIG. 5: Network formed by the BEDT-TTF molecules in 
the be plane of k-(BEDT-TTF)2Cu2(CN)3. Each grey ball 
represents the center of gravity of a BEDT-TTF molecule, 
with grey bonds indicating the intermolecular distances that 
were taken into account for the tight binding fit. BEDT-TTF 
dimers are highlighted by ellipses, and the paths forming the 
triangular lattice paths are shown with bold lines. The unit 
cell is marked as a rectangle. 



strength U/t estimated in this way monotonously falls 
by 7% as the temperature is increased from T = 5 K to 
T ^ 300 K. 

In order to rationalize the observed temperature de- 
pendence of the Hamiltonian parameters in ^ we ana- 
lyze the crystal structure in more detail. For that pur- 
pose, we determined the orientation of the BEDT-TTF 
molecules in space by finding the plane of the TTF part 
of the molecule and measuring its angle with respect to 
the be plane, cf. inset of Fig. [T] This yields the inclina- 
tion of the BEDT-TTF molecules against the anion layer 
shown as squares in Fig. [7] and the intradimer distance 
shown as circles (see also Table|I|. Both quantities show a 
nearly monotonous increase over the studied temperature 
range. The decreasing intradimer distance explains the 
increase of the intradimer hopping integral ti with de- 
creasing temperature. Apparently, the nonmonotonous 
evolution of the overlap integrals t, t', especially the dis- 
tinct extrema in both quantities around 150 K to 200 K, 
has to be of different origin. 

The overall trend can be understood by considering 
the temperature dependence of the lattice parameters; 
from Fig. [5] it is clear that changes in the b lattice pa- 
rameter should have an impact on t' , while changes in 
c should affect t. In Fig. [7] (b), we see that the c/b lat- 
tice parameter ratio decreases with temperature down to 
T = 150 K, then increases with falling temperature until 
T = 20 K. This has an immediate impact on the degree 
of frustration t' /t as it should be approximately propor- 
tional to the c/b ratio. Indeed, comparison of Fig. [6] (b) 
and Fig.[7](b) confirms this expectation and thus explains 
the nonmonotonous temperature dependence of the frus- 
tation. 

An increasing frustration upon cooling, reaching t' /t 
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FIG. 6: Temperature dependence of Hamiltonian parameters 
of K-(BEDT-TTF)2Cu2(CN)3. 



values at low temperatures in excess of those at room 
temperature, is an interesting finding which may help to 
better understand the intriguing low-temperature mag- 
netic and dielectric properties of this material. Here we 
mention the distinct types of charge ordering, accompa- 
nied by dielectric anomalies, proposed in Ref. [H for the 
present material as a result of an increasing degree of 
frustration. 

The nonmonotonous evolution of the in-plane distor- 
tion c/b, which adopts a minimum around 150 K, might 
be related to the ordering of the ethylene groups, un- 
covered in the present study^i According to our struc- 
tural analysis, the fraction of ethylene groups, ordered in 
the staggered conformation, gradually grows from 77% at 
room temperature over 86% at 250 K to 93% at 200 K. 
For temperatures below 150 K, the ordering is complete 
within the accuracy of our analysis/refinement. The large 
step size of 50 Kelvin employed in this study does not 
allow to determine the ordering temperature more pre- 
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FIG. 7: (a) Orientation i9 and dimerization d of BEDT-TTF 
molecules as function of temperature. The inset illustrates 
how i9 and d are measured, (b) Ratio c/b between c and b 
lattice vectors as function of temperature. 

cisely. Likewise, we cannot say whether or not the order- 
ing occurs continuously or sets in abruptly. An argument 
in favour of the latter possibility might be derived from 
a small step-like feature revealed in thermal expansion 
measurements around 150 K, see the out-of plane data 
shown in Fig. 1 of Ref. fol. We stress, that a fully ordered 
staggered ethylene conformation below 150 K and a pro- 
gressively disordered state above 200 K, is fully consistent 
with the anomalous behavior revealed by ^H-NMR mea- 
surements between 150-200 Ki^ The strong increase in 
(TiT)^^ above 200 K was attributed to thermally acti- 
vated vibrations of ethylene groups,- We suggest that the 
nonmonotonous temperature dependence in t' /t might 



also be related to the drastic change in the thermopower 
around 150 KJ The thermopower is related to the energy 
derivative of the density of states at the Fermi leveli^ 
However, for a strongly correlated system like k-(BEDT- 
TTF)2Cu2(CN)3, density functional theory is not suffi- 
cient for the calculation of this quantity and more elab- 
orate many-body calculations - which are beyond the 
scope of the present study - have to be done. 

In summary, we performed an analysis of the tempera- 
ture dependence of the structural and electronic proper- 
ties of K-(BEDT-TTF)2Cu2(CN)3 by considering a com- 
bination of X-ray diffraction at various temperatures and 
density functional calculations. Our study shows that 
the temperature dependence of the structural parame- 
ters has significant infiuence on the electronic properties 
and results in a nonmonotonous behavior of the degree of 
frustration. Of special relevance is the increase of frustra- 
tion at low temperatures in comparison to the behavior 
at room temperature. 
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Appendix A: Additional crystallographic material 

Crystallographic data for the k-(BEDT- 
TTF)2Cu2(CN)3 structure at 5, 10, 100, 150, 200, 
250 and 300 K has been deposited with the Cambridge 
Crystallographic Data Centre as supplementary publica- 
tion nos. CCDC 850022 to 850028. Copies of the data 
can be obtained free of charge on application to CCDC, 
12 Union Road, Cambridge CB2 lEZ, UK (fax: (44) 
1223 336-033; e-mail: data_request@ccdc.ccmi.ac.uk). 

Fig. [5] illustrates geometry of the anion laywer and the 
slight deviation from planarity. Fig. |9] and Table |ll] con- 
tain lattice parameters determined at 10 K intervals be- 
tween 290 K and 100 K. 
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T(K) a(A) b(A) c(A) PC) V{k') 




16.248 8.655 13.516 113.56 1742 

16.231 8.650 13.491 113.62 1735 

16.235 8.647 13.483 113.79 1732 

16.226 8.652 13.475 113.80 1731 
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FIG. 8: The Cul coordination sphere at 5 K. The Cul is 

0.072 A out of the plane defined by N12, Cll and C12. Ther- TABLE II: Unit ceU parameters as a function of temperature 
mal elhpsoids are drawn at the 50% probabihty level. in the temperature range from 290 K down to 100 K. 
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FIG. 9: Plot of the unit cell parameters as a function of tem- 
perature between 100 and 290 K. 
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